IJSET 



International Journal of Scientific Engineering and Technology 
Volume No.4 Issue No.6, pp: 370-377 



ISSN : 2277-1581 
01 June 2015 



A Novel Dynamic Model for a Crane System with a Flexible Cable and Large 

Swing Angle 

1 *Mohammad H. Fatehi, 2 Mohammad Eghtesad, 3 Roya Amjadifard 

department of Electrical Engineering, Science and Research Branch, Islamic Azad 
University, Tehran, Iran. 2 School of Mechanical Engineering, Shiraz University, Shiraz, 
Iran. Engineering Faculty, Kharazmi University, Tehran, Iran. 
Corresponding author: Mh_f atehi @ srbiau.ac.ir 



Abstract:In this paper a novel dynamical model will be 
developed for an overhead crane with transverse vibrations 
of flexible cable while considering large angle of the 
cable/payload swing. The major problem addressed is 
expressing the dynamic interaction between the cables and 
other elements in the crane system which demands 
investigation into boundary induced cable vibrations and 
large swing angle. The developed model is achieved as an 
ODE model and it is appropriate enough to design suitable 
control systems. Thus, to perform the modeling, Rayleigh- 
Ritz discretization method is used to achieve an ODE model 
for transverse deflection of the cable having finite 
generalized degrees of freedom. Using the Euler-Lagrange 
formulation, the dynamical model of the crane system with 
a flexible cable is obtained. The simulations are performed 
based on OpenSees model and the developed model to 
validate the obtained model. To demonstrate the dynamic 
behavior of crane system, an illustrative example of the 
crane system with flexible cable moving a lightweight and 
heavyweight payload is investigated under various 
conditions. 

Keywords — Overhead crane system, Flexible cable, Cable 
vibration, Large swing angle, Transverse vibration 

I. Introduction 

The crane systems are extensively employed in a variety of 
applications in industries such as land and onshore/offshore 
construction sites, transportation industry, etc. The most 
common operation of a crane is point-to-point moving of a 
suspended load horizontally, using the cables and a support 
mechanism. The cables possess an inherent flexibility and 
can only develop tension; they do not offer resistance to 
bending moments or compressive forces. Such natural 
features certainly cause deflection in transversal direction of 
the cable and payload swing in crane systems. The suspended 
load in crane systems is always subject to swings happen by 
unskilled operators or by disturbances typically induced by 
motor drive transients, wind, and collision with objects so 
that it can cause lengthy transportation activities; even the 
swings may become large and reduce the safety. 
Abdel-Rahman and Nayfeh presented a detailed review of the 
challenges in modeling and control of the crane systems [1]. 
In most dynamical models, the effects of flexibility and 
weight of the suspended cable have been disregarded and the 



cable has been considered as a mass-less rigid-link or as a rigid- 
link including a point-mass [2-14]. In these studies, the payload 
swing has been assumed as the major dynamic motion of crane 
systems. Although such assumptions are usual, they are not 
genuine for many applications. In certain cases, especially 
when payload is lightweight and more importantly when cable 
is long, the effect of flexibility has to be taken into account. In 
these cases, the tension force is more dependent on the cable 
weight and also it will be varying along the cable. Under such 
conditions, the tension force, especially at the end of the cable 
is low so that the cable weight can possibly have more 
dynamical effects and transverse vibrations of the cable may 
administrate the behavior of the crane [15]. Thus, to utilize the 
crane systems in a particular application, in addition to the 
payload swing, the cable vibration should be suppressed within 
a given period of time for safety issues; therefore, development 
of an effective suppression control system is essential. In order 
to achieve these objectives, a more practical and accurate 
model with more details including interactive dynamics of the 
cable and the payload is required. 

A few studies have addressed the effects of cable flexibility in 
crane systems [16-26]. They have offered planar models for 
overhead cranes and assumed that the cable is perfectly flexible 
and inextensible. Also, they have assumed that the cable slop is 
small along the cable and also disregarded the swing angle; in 
other words, in these studies, in a crane system, the major 
dynamic motion of cable has been considered to be the cable 
vibrations. However, this is valid only for slow movements of 
the trolley or support mechanism and near the end of the 
traveling. These models might not be accurate enough when 
considering certain unavoidable physical and environmental 
conditions. The angular rotation of cable swing in most 
applications especially when a high speed traveling is required 
becomes large. This issue has been disregarded in nearly all of 
the previous studies on the crane system with a flexible cable 
[17-26]. In an infrequent study, an approximated model was 
establish by introducing the effect of small swing angle just in 
the trolley dynamics in a crane system with flexible cable [16]. 
Here, an attempt to extract an accurate dynamical model that 
exactly represents the dynamic interaction effects of the 
payload and suspended cable and no particular assumption is 
made regarding the major dynamics of a crane system and both 
cable vibrations and payload/cable swing will be considered in 
dynamic modeling of the crane system. Moreover, the swing 
angle of the cable/payload is not being limited to small angles. 
Thus, a more accurate dynamical model including both the 
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transverse vibrations and large swing angle will be obtained. 
To do so, an overhead crane system with a flexible cable is 
considered and during the crane travel, it is assumed that the 
cable is inextensible with no load hoisting or lowering. The 
cables can be assumed as one-dimensional continuum 
elements which can be modeled approximately, assuming to 
have only a few modes of vibration [27-29]. One frequently 
used method for approximate modeling of continuum 
elements is Rayleigh-Ritz discretization. This method 
produces an approximate solution in the form of a finite 
series consisting of known shape functions or trial functions 
multiplied by undetermined time-dependent "generalized 
coordinates" [28]. Accordingly, the cable transverse 
deflection can be modeled by a finite number of generalized 
degrees of freedom and a set of ODEs using Euler-Lagrange 
formulation, which permits the determination of the 
generalized coordinates. 

The purpose of this study is to develop a novel dynamic 
model for a crane system with a flexible cable while 
considering the cable transverse vibrations and large angle of 
the cable/payload swing. The major problem addressed is 
expressing the dynamic interaction between the cables and 
other elements in the crane system which demands 
investigation into boundary induced cable vibrations and 
large swing angle. To attain an ODE model, Rayleigh-Ritz 
discretization method is used for transverse deflection of the 
cable having finite generalized degrees of freedom. Using the 
Euler-Lagrange formulation, the dynamical model of the 
crane system with a flexible cable is obtained. To 
demonstrate the effectiveness of the control system, the 
numerical simulations are performed using commercial 
software. The simulations are performed based on OpenSees 
model and developed model to validate the obtained model. 
To demonstratethe dynamic behavior of crane system, an 
illustrative example of the crane systems with flexible cable 
moving a lightweight and heavyweight payload are 
investigated under various conditions. 

This paper is organized as follows: The system description 
and kinematic analysis are introduced in Section 2. In Section 




Figure 1: Coordinate frames and the schematic of the 
overhead crane system. An infinitesimal element ds at 
general point p. 



3, the energy functions are introduced and the dynamic 
equations of motion are derived using Rayleigh-Ritz 
discretization method and Euler-Lagrange formulation. The 
several numerical simulations are performed Section 4, to 
illustrate and verify the behavior of a sample flexible cable 
crane system. Finally, in Section 5, the conclusions are drawn. 

II. System Description and Kinematic Analysis 

An overhead crane system is composed of a support mechanism 
as a trolley and a flexible cable tied to the suspended payload. 
The swing motion of the payload and transverse vibrations of 
cable can be described in a planar form using two coordinate 
frames i.e., XZand XZ, seeFigure 1. A general point p along the 
cable is defined to describe the position of a general point on 
the cable. The point p is located on a small element of length 
ds with the left end placed at point a and the right at point b. 
As shown, three kinds of motion can be considered, i.e. crane 
traveling, swing angle and transverse vibration which are 
described as x, 9, and v(z~t), respectively. It is assumed that 
every point on the cable has two degrees of freedom; one is the 
transverse deflection v(z~t) around the Z-axis and the other is 
the swing angle, 9, which is not assumed to be small in this 
study. Let x and F x be the trolley position and trolley driving 
force, respectively. The parameters M t ,m pl p and g are the 
total trolley mass, payload mass, mass per unit length of cable 
and gravitational acceleration, respectively. The cable is 
supposed to be inextensible and the transverse deflection is 
assumed to be small. The payload is considered as a point-mass 
and the motion of the trolley on the rail is assumed to be 
frictionless. According to Figure 1, the position vector of the 
point p in the base frame XZ can be described as: 



cos(i9) sin(i9) 
- sin(i9) cos(i9) 



_ (cos{d) (z, t) + sin(i9) z + x\ 
V - sinftf) (z, t) + cosftf) z ) 



(1) 



sin(i9) (z, t) + cos(i9)z 
where, x,z are, respectively, the trolley position and the 
component of the general point p position along the Z-axis. 
v(z, t)is the transverse deflection of the general point p and 
R ($) is the rotation matrix, in which d — 9 + e is the total 
angular displacement of point p. As shown inFigure 1, 9 and e 
are defined around Z-axis respectively are the payload swing 
angle and the cable slope at point p. It is assumed that the 
deflection of the cable is small, thus, the angle £ compared with 
the angle 9 is small, and therefore^ = 9. The velocity of the 
point p can be obtained as follows: 
— sin(0) cos(0) 1 /v\ 
-cos(0) -sm(9)\\z) 

os(9) sin(0)l /v t \ ,x\ 
sin(0) cos(0)J \z t ) KqJ 
where, v t and z t denote the partial derivative of v and z with 
respect to time, respectively. Since it is assumed that the cable 
is inextensible, all points along the cable have the same axial 
velocity, thus, z t = z = 0, and, given these assumptions, the 
approximated velocity of point p is calculated as follows: 
x — 9 sin(0) v + 9 cos(0) z + cos(0) v t N 
—9 cos(0) v — 9 sin(0) z — sin(0) v t 



+ 



(2) 



(3) 
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(4) 



(5) 



To attain an ODE model describing the transverse deflection 
of a cable with finite degrees of freedom (modes), the 
Rayleigh-Ritz discretization method can be used [28]. Based 
on Rayleigh-Ritz discretization, the spatial function v(z',t) is 
approximated as the finite sum of shape functions 0 ; (z) 
multiplied by the time-dependent generalized coordinates 
ji i (t) as follows: 

m 

V (z, 0 = £ Hj (00; (i) = P T (z)Kt) 
7=1 

where, P(z) and p(t)are: 
P(£) = [<Mf) 02 (z) - 0 m (z)] r 

Thus, the partial derivatives of v(z~t) with respect to time 
and zcan be calculated as: 
v t (z,t) = P T p 

d d 2 (6) 

dz dz A 
To choose the shape functions 0 ; (z), the boundary 

conditions must be satisfied. It is assumed that the boundary 
conditions are v(0, t) = v{C,t) — 0 so the shape functions 

can be chosen as 0, (z) — sin (y^) ' V = 1'2<3, To 
obtain the velocity of point p, the following relationship can 
be used. 

(± - 0 sin(0) P T p + 0 cos(0) z + cos(6l) P T p" 



Tp \ -B cos(6) P T p - G sm(G) z - sm(B) P T p 1 
In the following subsections, these formulations will be used 
to obtain the potential and kinetic energies required when 
extracting the dynamic equation based on Euler-Lagrange 
formulation. 



III. Dynamic Equations of Motion 

Let an infinitesimal element ds be on the point p along the 
cable, as shown inFigure 1, where, ds — V(dz) 2 + (dv) 2 . It 
is assumed that the cable longitudinal elastic deformation is 
neglected. The cable transverse deflection is also assumed to 
be small, thus, ds = dz. So, the kinetic energy of the crane 
system can be expressed as: 



(8) 



K= \p\ f p T r p dz + X -M t x 2 + fn p {r p T r p ) 
Jo 

where, the first term denotes the kinetic energy of the cable 
and the other two terms are related to the trolley and the 
payload, respectively. The potential energy of the system is 
due to the payload and the cable gravity effects, and the cable 
transverse deflection. The potential energy due to the cable 
transverse deflection can be obtained as: 



U = pg(z - 
Jo 



zcos (0) + v sin(60). dz + i T{z)wldz 



r 

Jo 



+m p g(£ - {cos(.6) + v{£, t) sin(0)) 



where, , the first term is the potential energy due to the cable 
transverse deflection, second and third terms is potential 
energy due to gravity effects of the cable and payload, 
respectively. 7(z)is the tension force along the cable, w$ is 
the partial derivative of v[z, t) with respect to z, and I is the 
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cable length. Since the cable is hung under its own weight and 
also weight of the payload, the tension in the cable at point p 
can be approximately obtained as: 

T(z) = (m p g + pg(l - z)) cos(0) (10) 

It has been assumed that both the deflection and slop of the 
cable aroundz-axis are small, thus, cos(e) = 1, and sin(e) = 
£ = tan(e) = v f , so, 

cos(i9) = cos(0 + s) = cos(0) — tan(e) sin(0) n 
= cos(0) - i; z -sin(60 ( } 

Thus, the cable tension can be obtained approximately as: 

T(z) = T 0 (z) cos(0) - r 0 (z) sin(0) v z - (12) 

where, T 0 (z) — (m p + pi — pz)g. Therefore, the total 
potential energy can be approximately written as following 
where, U 0 = \pgt 2 + m p g{ is a constant value. 

U = - r 0 (z) cos(#) v\dz + pg{z — z. cos(#)). dz 
Jo Jo 



(13) 



+m p gi — m v glcos{6') 

= -\ T 0 (z) cos(#) v\dz — pg I zcos(9).dz 
2 Jo Jo 
—m p g£cos(8) + U 0 

To obtain crane's equations of motion based on the Euler- 
Lagrange formulation, the system Lagrangian as a function of 
the trolley motion, cable swing (pendulum) motion and cable 
(discretized) deflection will be used. This Lagrangian function 
is given as: 
L — K — U 

\p( f p Tf p) ~ K m p9 + pgi)cos(8)v] 



'P 'Pj 
+ 1 



(l + it7f)p 5 zcos(0) 
+\M t x 2 + \m p {r p T r p ) f=[ + m p< gtcos(0) 

1 T 



! dz 



Un 



(14) 



2 1 2 

Assuming that, q — (x, 9, p) ' £ R n in which p — 
{p 1 ,p 2 , ■•■ ,p m )eR m is generalized coordinates of cable, and 
n — m + 2. To attain the dynamic equation, Euler-Lagrange 
equation is used as: 

d (dL\ dL 

dt \dq) dq 

By inserting the Lagrangian function (14), L(q,q), into Euler- 
Lagrange equation (15), the dynamic equations of the crane 
system will be obtained in a matrix form as: 

M(q)ij + C(q, q)q + G{q) = Wu (16) 

where, W — [1 0 ••• 0] r , u = F x is the trolley driving 
force and q — (x,9 ,p) T £ R n in which 
p — (pi,P2, "• ,p m ) e P m is generalized coordinates of cable, 
and n — m + 2. Matrix M(q) is the total inertial matrix, 
C{q,q)q represents the Coriolis and centripetal forces and 

(9) G{q) — -j- is the potential force due to stiffness and gravitation 

effects in the system. The symmetric inertia matrix, M(q), can 
be assembled and simplified as: 

a 0 cii cos(0) cosi;p?)a 2 



M(q) = 



a-t cos(0) a 4 + p T a 3 p 
cos(0) a\ al 



(17) 
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C(q,q) = 



G(q) 



(18) 



Also, matrix C{q, q) and vector G{q) can be assembled as: 
0 c 12 — a 2 sin(0) 8 
0 n T a 3 ii 9n T a 3 
0 -0a 3 /i 0 
( 0 
g 1 sin(0) + sm(6)ii T g 2 ii 
\ cos(6)g 2 ii 
in which, c 12 = — % sin(0) 0 — a 2 fi cos(0) 0 — sin(0) a 2 fi 
and the respective parameters are: 



a n = 



a 2 



M t + m p + pgi 

t 

P T dz 

o 



Jn 



a 3 



a 5 



= p \ (P T P)dz 
Jo 

= p I zP T dz 
Jo 



9i = {m v i + ip^ 2 ) g 



(19) 



52 



■* n 



These nonlinear coupled ordinary differential equations 
describe the dynamic motions of the whole crane system with 
flexible cable. It is clear that q is the configuration variable 
vector of the system with q 1 — x is the actuated variable and 
q 2 — (0, ji) T is the un-actuated variable vector of the system. 
Since there are three configuration variables to be controlled 
with only one actuated configuration variable, the flexible 
cable crane system is an underactuated system. There are 
some general properties of inertia matrix M{q), note that this 
matrix is symmetric, and positive definite for all q. 
Calculating M — 2C is a skew-symmetric matrix which has 
an important property. 

IV. Numerical Simulation to Verify the Model 

In this part, several numerical simulations are performed to 
illustrate and verify the behavior of a sample flexible cable 
crane system with a lightweight and heavyweight payload. 
To do so, an OpenSees computer model is used to compare 
and validate the developed model. Since the crane system's 
degrees of freedom are cable deflection, payload/cable swing 
angle and trolley movement, time histories of these motions 
will be plotted to illustrate the behavior of the sample crane 
system. To collect the simulation results, the obtained ODE 
model with four finite modes, as described in the previous 
section for a cable with suspended payload, is used. Also, a 
model with 100 dimensional finite elements which is built in 
OpenSees software is applied to investigate and compare the 
developed model. Assuming that there is no load hoisting or 
lowering during crane travel, two cases are simulated to 
verify the obtained model. First, a crane with a flexible cable 
and lightweight payload and, then, the same crane with a 
heavy payload are examined. For example, the trolley mass is 
10 (kg), cable length is 5 (m) and its mass in per unit length 



is 0.62 (kg/m), the mass of different payload in two cases are 
chosen 2 (kg) (light payload case) and 50 (kg) (heavy payload 
case). It is assumed that the crane system is initially at rest. To 
excite the cable-payload assembly, the same acceleration signal 
as a positive and negative pulse is imposed on the trolley 
(Figure2). To show traveling displacement and transverse 
deflection of the cable, five chosen generalized points with the 
same distance along the cable are considered, where traveling 
displacement related to the end of cable is payload motion. 
Assuming that the crane system is in the rest condition before 
moving trolley, the simulation results are plotted in Figure and 
Figure, to illustrate horizontal displacement of the cable and 
payload in line with the X-axis (horizon) for two case cranes. 
As seen, the dynamic behaviors of crane system based on 
developed model and OpenSees model are close together and 
the obtained model, according to the previous section with four 
finite degrees of freedom yields good results. To investigate the 
other dynamic behavior of crane system, another simulation is 
performed. In another way, the same crane system is excited 
via the same trolley force as a positive and negative pulse 
(Figure3). The dynamic behavior of crane system is shown in 
Figure4 and Figure7 to Figure9 for two cases with lightweight 
and heavy payloads. To determine the transverse deflection of 
the cable, the first four modes of vibrations are considered. The 
time histories of the amplitudes of these first four modes are 
plotted in Figure4 and Figure6. This figure clearly shows that 
the amplitudes are diminishing with the number of the mode, 
and thus, justifies using the first four modes for the 
approximated system in our simulations. In Figure8 and 
Figure9 the horizontal displacements of the trolley and the 
payload along X-axis is displayed. Also, Figure presents the 
payload swing angle. The results show that the swing angle can 
become large and it cannot be disregarded. To investigate the 
effectiveness of the payload mass, the time histories of 
transverse deflections of the cable is shown in 
Figure5andFigure7. The results both illustrate the crane system 
behavior and the amount of effectiveness of payload weight on 
cable vibrations. As shown, when the payload is lightweight, 
the transvers deflection of cable and its effectiveness on 
payload swing become more because the tension force along 
the cable is not large enough, especially at the end of the cable 
which is allowed to have more open dynamical behavior. 
Obviously, when a heavy payload is tied to the end of cable, 
tension force on the entire length of the cable is large, so the 
transvers deflections of cable are small and the effectiveness of 
cable dynamic on payload can be ignored. Under such 
conditions, cable can be assumed as a rigid-link. 

V. Discussion and Conclusion 

This study has developed a comprehensive and convenient 
dynamical model for an overhead crane with a flexible cable 
while allowing large payload/cable swing angles. The cable 
was assumed to be a one-dimensional continuum element. To 
discretize the governing equations of this continuum, the cable 
transverse vibrations were approximated with finite generalized 
degrees of freedom based on Rayleigh-Ritz discretization 
method and a nonlinear model for the whole crane system was 
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developed using Euler-Lagrange equation. As it can be seen, 
the dynamic responses of the crane system based on the 
developed model and OpenSees model are very similar and 
the obtained model generates decent results. The results both 
illustrate the crane system behavior and the amount of 
effectiveness of payload weight on swing angle and cable 
vibrations. As shown, when thepayload is lightweight, the 
transvers deflection of cable and its effectiveness on payload 
swing become more significant. Thus, considering both 
dynamic behavior of the swing angle and transverse 
vibrationsare important and cannot be disregarded. But, when 
a heavy payload is tied to the end of cable, the transvers 
deflections of cable are small and the effectiveness of cable 
dynamic on payload can be ignored, thus, under such 
conditions, the major dynamic behavior of cable is swing 
angle and cable can be assumed as arigid-link. 
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Figure5:The comparison of the simulation results using the 
OpenSees model and developed model. The horizontal 
displacements of five points, located on the cable at lm, 2m, 
3m, 4m and 5m, around vertical axis for a crane system with 
2kg payload 



IJSET@2015 



Page 374 



IJSi 



International Journal of Scientific Engineering and Technology 
Volume No.4 Issue No.6, pp: 370-377 


























































































Opensees model result 

Developed ODE model result 

















Time (sec) 



































































































Op 


?nsees model re 


ult 










Developed ODE model result 










1 1 





Time (sec) 

































































Opensees model r 


suit 

del result 




Developed ODE mc 















Time (sec) 

Figure6:The comparison of the simulation results using the 
OpenSees model and developed model. The horizontal 
displacements of five points, located on the cable at lm, 2m, 
3m, 4m and 5m, around vertical axis for a crane system with 
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Figure4: The time histories of the amplitudes of the first four 
vibration modes of the same crane system with a2kg payload 
under the exciting input force 
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Figure5: Time histories of transverse deflections of the four 
selected points along the cable for the same crane systemwith 
a2kg payload under the exciting force 
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Figure7:Time histories of transverse deflections of the four 
selected points along the cable for the same crane systemwith 
a50kg payload under the exciting input force 
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Figure8: Time histories of horizontal displacements of the 
trolley and the payload for the same crane system with a2kg 
payload under the exciting input force 
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Figure9: Time histories of horizontal displacements of the 
trolley and the payload for the same crane system with a50kg 
payload under the exciting input force 
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